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Abstract 

Using the collinear QCD factorization approach, we study the single-transverse-spin dependent 
cross section Aa(S±) for the hadronic production of two jets of momenta P\ = P + q/2 and P2 = 
— P + q/2. We consider the kinematic region where the transverse components of the momentum 
vectors satisfy P± 3> q± 3> Aqcd- For the case of initial-state gluon radiation, we show that at the 
leading power in q±/P± and at the lowest non-trivial perturbative order, the dependence of Aa(S±) 
on q± decouples from that on P± , so that the cross section can be factorized into a hard part that 
is a function only of the single scale P±, and into perturbatively generated transverse-momentum 
dependent (TMD) parton distributions with transverse momenta k± = 0(q±). 
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1. Introduction. Single-transverse-spin asymmetries (SSAs) in high-energy hadronic 
reactions with one transversely polarized hadron were first observed more than three decades 
ago 1]. The SSA is defined as A N = (a(S±) — <y(—S±))/(a(S±) + a(—S±)), the ratio of the 
difference and the sum of (differential) cross sections when the hadron's spin vector, S±, is 
dipped. Recent experimental measurements of SSAs both in polarized hadronic collisions 
2I, y] and in semi-inclusive lepton-nucleon deep inelastic scattering (SID IS) J4] have renewed 
the interest in investigating the origin of SSAs in Quantum Chromodynamics (QCD) jlj]. 

It is believed that some SSAs are a consequence of the partons' transverse motion inside 
the polarized hadron. The momentum scale of this transverse motion is a typical hadronic 
scale, (k±) ~ 1/fm ~ Aq CD . For observables with only one hard scale Q 3> Aq CD , the SSA 
should be proportional to (k±)/Q . Such observables only probe an averaged effect of 
the partons' transverse motion. However, for observables characterized by more than one 
physical scale, SSAs may directly probe the partons' transverse motion. For example, in 
the case of Drell-Yan hadronic production of a lepton pair of large invariant mass Q and 
transverse momentum q± <C Q, the pair probes the (anti-) quark's transverse motion at the 
scale q±, while the invariant mass Q of the pair sets the hard scale of the collision In 
this letter, we study the SSA in hadronic production of two jets: A(Pa,S±) + B(Pr ) — > 
Ji(Pi) + J 2 (P 2 ) + X, with the jet momenta Pi = P + q/2 and P 2 = -P + q/2 0, Q,kl|. 
Unlike in the Drell-Yan process or in SIDIS, the SSA in di-jet production can be generated 
by both initial- and final-state interactions. We emphasize that measurements of the SSA 
for di-jet production have begun at RHIC 12], complementing the measurements in SIDIS. 

We are interested here in deriving a QCD formalism for the single transverse-spin de- 
pendent cross section, Aa(S±) = (a(S±) — cr(—S±))/2, that is valid in the kinematic region 
P± ^ Q± ^ Aqcd, where P± and q± are the transverse components of the momenta P and q, 
respectively, so that the SSA provides direct information on the partons' transverse motion. 
We first consider the region P± 3> q± 3> Aqcd, where both observed momentum scales are 
much larger than the typical hadronic scale Aqcd- We calculate Aa(S±) in terms of the 



collinear QCD factorization approach, which is expected to be valid in this region 131 ] . In 
this approach, the incoming partons are approximated to be collinear to the corresponding 
initial hadrons, and the leading-order partonic processes produce two back-to-back jets with 
zero momentum imbalance. The di-jet momentum imbalance, q± = P\± + P 2 j_, has to be 
perturbatively generated by radiating an additional hard parton. 
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(a) (b) (c) 

FIG. 1: Sample diagrams for quark-quark scattering contributing to Aa(S±) through an initial-state 
interaction (a), and through final- state interactions with jet Pi (b) and jet P2 (c). 

In this letter, we concentrate on the physics issues related to di-jet production, and re- 
strict ourselves to the case that the leading contribution in the expansion of the partonic 
scattering in q±/P± involves a hard qq' — > qq' subprocess. We consider the jet imbalance 
q± to be generated by gluon radiation off the initial quarks. We derive the corresponding 
leading-order contribution to Aa(S±), and demonstrate that the perturbatively calculated 
partonic parts can be further factorized into a single-scale (P±) hard part and perturbatively 
generated transverse-momentum dependent (TMD) parton distributions with transverse mo- 
menta k± = 0(q±). We find that the complete contributions from all other partonic sub- 
processes at the leading order have the same factorization property. These will be discussed 
in a forthcoming publication [14J]. The factorization of the physics at scale P± from that at 
scale q± that we find when q± <C P±, is consistent with a more general TMD factorization 
formula for the SSA in the di-jet momentum imbalance. 

2. Single transverse-spin dependent cross section. When both P± and q± are much 
larger than Aqcd, a nonvanishing single transverse-spin dependent cross section Aa(S±) is 
generated by the Efremov-Teryaev-Qiu-Sterman (ETQS) mechanism in the collinear 

factorization approach. The calculation of Aa(S \ ) then requires to evaluate partonic pro- 
cesses with 3-parton initial- and final-states [13|. In Fig. [TJ we show generic diagrams for 
the quark-quark scattering channel that contribute to Acr(S±) through initial- and final- 



state interactions with the gluon o 
required for a nonvanishing SSA 



momentum k g , which is needed for generating the phase 



a a. 

mentum k' into 

the final state generates the jet imbalance q±. The blob in the center represents tree-level 
Feynman diagrams with the given initial- and final-state partons. In the ETQS formalism, 
the contribution of the subprocess (g)qq' — »■ qq'g to Aa(S±), shown in Fig. [TJ is generically 
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given by 



dAa(S±)( qq >) 



I 



dx' 



— dx\dx2 Tp(xi, X2) q'{x') 



dyidy 2 dP]_d 2 q± 



X' 



1 



(1) 



x 



16s(2vr) 4 



where y\ and y 2 are the rapidities of the two jets, s = (Pa + Pb) 2 , and Tt represents a 
partonic hard part. q'(x') is the usual quark distribution at momentum fraction x' in the 
incoming hadron B. x\ and X2 are the momentum fractions of the quarks from the polarized 
hadron A on the two sides of the cut shown in Fig. [TJ and Tp(x\,X2) is the corresponding 
twist-three quark-gluon correlation function, extracted from the lower blob in the figure 



where C is the proper gauge link to make the matrix element gauge invariant, and where 
the sums over color and spin indices are implicit. In Eq. and the rest of this paper, the 
dependence on factorization and renormalization scales is suppressed. 

Equation <Q applies when q± ~ P±. Our goal is now to investigate the leading structure 
that emerges from Eq. ([I]) when qj_ <C P±. In this limit the gluon of momentum k' is 
radiated either nearly collinearly from one of the external quark legs and/or is soft. In the 
present work, we only discuss collinear emission by one of the initial quarks. This radiation 
is the most interesting from the point of view of studying the factorization properties of 
the cross section at small q±, because TMD factorization can only hold if the initial-state 
collinear radiation leads to a certain specific structure, as we shall discuss below. Since 
we are considering the production of jets (as opposed to that of two specific hadrons), 
collinear radiation from final-state quarks becomes part of the jet and will not produce 
leading behavior in q±/P±. On the other hand, large- angle soft gluons produced by the 
interference of initial- and final state radiation may give leading contributions, through a 
so-called soft factor. We leave the detailed study of the soft factor to future work, but will 
briefly return to it later. 

As we mentioned above, the strong interaction phase necessary for a nonvanishing Aa(S±) 
arises from the interference between the imaginary part of the partonic scattering amplitude 
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FIG. 2: Sample diagrams for the soft-pole (1-3) and hard-pole (4-6) contributions. We have indi- 
cated the momentum k' of the radiated gluon that produces the jet imbalance q±, and the momentum 
kg of the additional gluon from the polarized proton. The pole of each diagram is taken from the 
propagator with a short bar. Diagram (I) is an example of an initial-state interaction, and (2) and 
(3) show final-state interactions. Diagrams (4-6) show the complete set of diagrams for one hard 
pole, indicated again by the short bar. 



with the extra polarized gluon of momentum k g = x 9 Pa and the real scattering amplitude 
without the gluon in Fig. [TJ The imaginary part comes from taking the pole of the parton 
propagator associated with the integration over the gluon momentum fraction x g — Xi — X\. 
For a process with two physical scales, P± and q±, tree scattering diagrams in Fig.[T]have two 
types of poles, corresponding to x g = ("soft-pole") [j| and x g ^ ("hard-pole") ^]. With 
the extra initial-state gluon attachment of momentum k', there are many more diagrams 
that contribute to Aa(S±) in comparison to the spin-averaged cross section In Fig. [2l 
we show some sample diagrams that give the leading soft-pole (1-3) and hard pole (4-6) 
contributions to Aa(S±) at q± <^ P±. By using the power counting technique [la ], we are 
able to classify all Feynman diagrams into different groups 14J. For example, diagrams (1-3) 
only make leading contributions when the momentum k' of the radiated gluon is parallel to 
Pa, whereas diagrams (4-6) can give a leading contribution when k' is either parallel to Pa or 



5 



to Pb- Beyond that, the calculation of the soft-pole and hard-pole contributions to Aa(S±) 
follows the same procedure as introduced in Ref. 8| for the SSA in Drell-Yan and SIDIS, 
because the kinematic limit considered is similar. In order to extract the contributions to 
A<t(iSj_) in Eq. (jTj), we need to convert the extra gluon field operator in the hadronic matrix 
element of the polarized hadron in Fig. [T] to a field strength operator in the definition of 
Tf(x\,X2) in Eq. (JSJ). Working in Feynman gauge, we first give the initial-state collinear 
partons from the polarized hadron a small transverse momentum, fcj = X{Pa + ki± with 
i = 1,2, and then expand the calculated partonic scattering amplitudes around ki± = 0, or 
equivalently, k g ± = k2± — k\± = 0. The contribution to Aa(S±) arises from terms linear in 
k g ±. After summing up all contributions [ijj], we obtain the total leading-power contribution 
to Aa(S±) from the (g)qq' — > qq'g partonic subprocess in the q±/P± expansion: 

dAa(S_L)( qq f) 



dy 1 dy 2 dP]_d 2 q_ l 
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dx 



X 



(<?1) 2 

A + T F (x a ,x a ) 



a/3 QOL n P 

J ± i l± rrSivers 

qq'->qi' 



dx' 



x 



B 



(3) 



where 



A 



9^, , 1 

x df FM W r 



+Tp(x, x) 



1 



2N r 



B 



+T F (x,£x) 



q'(x') C F 



2N r 



[i+e] 



(4) 



1+j 
1-^ 



z /2 



where ^ = x a /x and = Xb/x' with x a = ^ (e Vl + e V2 ) and Xb = ^ (e~ Vl + e~ 
above relations are regularized at the integration limits by "plus" -distributions 
single-scale partonic hard part is given by 
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H^ s qql (s,t,u) 



a 2 s 7i 



AN? 



2(s 2 + u 2 

i 2 



(5) 

The 
The 



(6) 



where the use of the superscript "Sivers" will become clear in the next section, and 
where the partonic Mandelstam variables are given as s = x a XbS, t = —P 2 (e y2 ~ yi + 1), 
and u = —P]_(e yi ~ V2 + 1). We note that to the leading power in q±/P±, we have 
P 1 = P ± (e yi / V2,e- yi /V2,l) and P 2 = P± (e^/v 7 ^, e' m j\f2, -1) for the jet momenta in 
light-cone coordinates. 



The hard part for the partonic channel we have considered, ,, is very similar to the 



spin-averaged partonic differential cross section da/dt 14j. The only difference is the color 



factor in square brackets. In fact, the color factor for H^^ B , is equal to a sum of three 
color factors: C / + C 'f 1 + C p 2 , corresponding to color factors of scattering amplitudes when 
the initial-state gluon of momentum k g is attached to the initial-state incoming quark of 
momentum the final-state quark with P-u or the final-state quark with P2, respectively. 
For the (g)qq' — > qq'g subprocess, we have 14 1 Cj — — ^hr, Cp 1 = — jhr, and Cp 2 — A ~~ 



2iV£' ^-fl 4/Vg,' " ^ ^2 4AT2 



As a result, the final-state interaction with the jet of momentum P 2 dominates, and the 
overall color factor C\ + Cp 1 + Cp 2 has a sign opposite to that of Cj alone. 

3. Factorization in terms of TMD distributions. When q± <C P±, the di-jet 
production and Drell-Yan process at low transverse momentum (qp Q) considered in 
Ref. [8| share very similar kinematics. The difference is that the Drell-Yan process has only 
initial-state interactions while for di-jet production both initial- and final-state interactions 



are present. It is known [171 . Il8| that when q± Q, the Drell-Yan cross section in leading 
order in q±/Q can be calculated from a generalized QCD factorization formula involving the 
TMD parton distributions. It is natural to ask if the generalized QCD factorization formula 
can be extended to the di-jet cross section for q± <C P±. 

When q± ^> Aqcd? TMD parton distributions can be calculated in perturbative QCD 
(pQCD) from hard radiation and parton splitting [8j. The unpolarized TMD quark distri- 
bution is well known and given by [8] 

where B is given in Eq. (jSJ), and where the ellipses denote terms proportional to 8(1 — £) and 
terms generated by gluon splitting. The transverse-spin dependent TMD quark distribution, 



qT{x a ,q±), known as the Sivers function 19], can also be calculated in pQCD and is given 



in terms of the twist-3 quark-gluon correlation function Tp as 

= [.4 +...], (8) 

where A is given in Eq. Mp is a hadron mass scale introduced to keep q(xb,q±) and 
qr{xa, q±) at the same dimension, and where the ellipses as above denote terms proportional 
to 5(1 — £) and terms not relevant to the following discussion. In Eq. (jSJ), the superscript 
"SIDIS" indicates the Sivers function for the SIDIS process, which has an opposite sign from 



that for the Drell-Yan process given in Ref. [8j, due to the difference in the directions of the 
gauge link that defines the TMD quark distributions 
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23(]. 



One of the important features of the di-jet cross section calculated above to leading order 
in q±/P± is the separation of two observed physical scales, P± and q±. We can use Eqs. (j7|) 
and dSJ) to rewrite the single transverse-spin dependent di-jet cross section in Eq. (j3J) in terms 
of a single-scale hard factor, H Slvers , which is a function of P±, and of the g_i_-dependent 
perturbatively generated TMD parton distributions: 

dAa(S ± ) {qql) 



dyidy 2 dP 2 d 2 q± 



a/3 na 8 rrSh 

e o ± q ± n , 



fx>9x>A QC D 



lvers 
qq'^qq' 
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M(x h )) {x a qf DlS (x a ,q ± )) 
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4" ( x bq'{x b ,q±))(x a T F (x a ,x a )) 
1± 



Using the leading order relation 23] 

1 



d 2 k ± k 2 ± qf Dls (x,k 



-T F (x,x) 



(10) 



we find that our result in Eq. (jUJ), when the contributions from all other partonic subprocesses 
at the same order [l^] are added, is consistent with the leading-order term of a more general 
factorization formula in terms of TMD parton distributions: 

e^Slq^ 

ab 

x a q^ IS (x a ,k 1± )x b fl lDls (x b ,k 2± ) 



dAa(S ± ) 
dy 1 dy 2 dP\d 2 q 1 

h± ■ q± 



[ d 2 k 1± d 2 k 2± d 2 Xj 



X ■ 
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X [Sab-*cd{^J 



(11) 



H!^(PT)] C 6®(ku. + h± + Xjl - q±] 



where apart from the functions already given, IDIS denotes the unpolarized TMD quark 
distribution and S a b^cd is the soft factor mentioned above Because of the color 

flow into the jets, the product of the soft and hard factors will involve a sum over separate 
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25], which has been represented by 



color amplitudes in the full factorization formalism 
a trace [ ] c in color space in the above equation. Our calculation of initial-state collinear 
gluon radiation described above would not be sensitive to this complexity of the color flow, 
but we emphasize that the definition of the parton distributions cannot be affected by it 



s 



In Eq. (fTTj) . we have chosen TMD parton distributions defined in SIDIS because of the 
dominance of final-state interactions. Choosing TMD parton distributions defined according 
to the Drell-Yan process would change the sign of the partonic hard factors, but not affect 
the overall sign of the physical cross section. Based on our explicit calculation here and the 

nn 

generalized factorization property of the Drell-Yan process at low q± |17l . 1181 ] , and because of 
the similarity in kinematics between two processes, we expect the generalized factorization 
formula in Eq. (fTl~j) to be valid for describing the single-transverse-spin dependent cross 
section for the di-jet momentum imbalance in hadronic collisions in the kinematic region 
where P± > q± > Aqcd- 

A key feature of the factorization is that the perturbatively calculated short-distance hard 
factors should not be sensitive to details of the factorized long distance physics. We tried, as 
a test, to derive all short-distance hard factors by using this factorization formula and the 
Brodsky- Hwang- Schmidt model for SSAs [20]. We were able to recover all hard-scattering 
factors in this way. We note that the same hard factorsi& r c s d 

as above were also found for 
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A detailed comparison between 



the weighted (integrated) SSA for di-jet production in 
the approach of lOj and ours will be presented in Ref . |14| • An all-order proof of the above 
factorization formula, if correct, is still needed and is beyond the scope of this paper. 

4. Summary. We have studied the single-transverse-spin dependent cross section 
Aa(S±) for di-jet production momentum imbalance in high-energy hadronic collisions in 
a kinematic region where P± 3> q± ^> Aqcd, and calculated contributions from both initial- 
and final-state interactions. At the leading order in q±/P±, the q± and P± dependences in 
our calculated results are decoupled and can be factorized into a single-scale hard factor that 
depends on P±, and into perturbatively generated TMD parton distributions. This factor- 
ization occurs for each partonic channel. Overall, final-state interactions turn out to give the 
dominant contribution to Aa(S±) [14]. We therefore expect the SSA in di-jet production to 
have the same sign as the Sivers asymmetry in SIDIS. 

We have found that our results are consistent with a more general TMD factorization 
formula, given in Eq. ( ITTi) . which we propose to be the correct approach to describing the 
single-transverse spin asymmetry in di-jet production at hadron colliders when P± ^> q±. 
We emphasize that obviously the result of our first-order calculation is not able to actually 
prove this factorization, but should rather be regarded as a "necessary condition" for such 
a factorization to hold. A full proof remains an important challenge for future work. Also, 
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we recall that we have limited ourselves to the case of collinear initial-state gluon radiation. 
The effects of large-angle soft-gluon emission have been neglected, even though we have 
indicated their likely role in Eq. fill I) . A proof of TMD factorization in this process would 
naturally incorporate a study of this soft factor. 

If the proposed factorization formula is valid, the di-jet momentum imbalance at RHIC 
could be described by the same TMD parton distributions as those used to describe the SIDIS 
and Drell-Yan processes. Because both initial- and final-state interactions are present, the 
di-jet momentum imbalance is sensitive to different short distance dynamics, and it will 
be an excellent process to test QCD factorization and the universality of the TMD parton 
distributions. In addition it should give valuable information on the partons' transverse 
motion in the nucleon. 
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